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0\ \ Abstract 

We study the vrvrA^ system by solving the Faddeev equations, for which the input two- 
body f-matrices are obtained by solving the Bethe-Salpeter equation in the coupled channel 
formalism. The potentials for the vrvr, vrA^ sub-systems and their coupled channels are obtained 
Cji ' from chiral Lagrangians, which have been earlier used to study resonances in these systems 

successfully. In this work, we find a resonance in the vrvrA^ system with a mass of 1704 — z375/2 
MeV and with quantum numbers I = 1/2, J'" = 1/2"'". We identify this state with the 
A^*(1710). This peak is found where the energies of the vrvr sub-system fall in the region of 
^ ■ the a resonance. We do not find evidence for the Roper resonance in our study indicating 

^D , a more complex structure for this resonance, nor for any state with total isospin I = 3/2 or 

^ : 5/2. 

^, 

O ■ 1 Introduction 

oo 

^^ ' The excited states of the nucleon have been studied extensively theoretically as well as ex- 

perimentally. This is evident from the fact that many of these states, especially those in the 

^ , energy region below 1750 MeV, have been assessed either three or four stars by the particle 

$H ' data group (PDG) [Ij . Even then, there are some resonances in this low energy region which 

still need unanimous agreement on their characteristics or existence, e.g., the J'^ = 1/2"'" 
resonances in the isospin 1/2 domain. The A^*(1440) or Roper resonance is a subject of 
continuous debate and the existence of the A^*(1710) is even questioned. The quark models 
face difficulties in reproducing both these states [21 [Sj [4]. In case of the A* (1710), some 
partial wave analyses [U [6] do not find any pole corresponding to it, while others claim a 
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clear manifestation of this resonance [7l [HI [U [101 • On the other hand, the authors of [llj claim 
an indisputable existence of the A^*(1710) from their study of the vrA^ -^ rjN reaction in the 
coupled channel formalism and suggest that the status of this resonance should be improved 
from three-star to four-star. 

Another controversy about the A^*(1710) started after the finding of a narrow peak in the 
'J A — > {K^n)X reaction at LEPS [12j, suggesting the existence of a pentaquark state which 
some groups associated to a SU(3) antidecuplet to which the A^*(1710) would also belong 
(see, for example, [I3l[l4j). In order to be compatible with the @~^, the A'"* (1710) is required 
to be narrow. However, the width of this resonance is not known precisely, with the widths 
listed in [T] ranging from ~ 90 - 480 MeV. The authors of [15] re-analyzed the vrA^ — > KA 
reaction and found that a narrow width of the A^*(1710) [1] was incompatible with the data 
and proposed the existence of another narrow resonance in this energy region. The partial 
wave analyses group who do not find a pole for the A^*(1710) suggested to look for other 
resonance in this energy region as a possible narrow, non-strange partner of the 9^ |161 [T7] . 
The debate on this issue has continued with new analyses which do not find a signal for the 
^+, as a consequence of which, the case for this state has weakend (see p!8] for a review). 

In case of the Roper resonance, which is the lowest excited state of the nucleon and, hence, 
in the simplest quark model should be expected to be a 3-quark state with a radial excitation 
of a quark, alternative descriptions, like a 3-quark-gluon structure [H], a quark core dressed 
by meson clouds [20], a dynamically generated resonance from interaction of mesons and a 
baryon [2l], etc., are posed in order to reproduce its properties. 

Looking at the characteristics of both these 1/2"'" resonances in [1], i.e., a large branching 
ratio for the vrvrTV decay channel, (~ 30-40 % for the iV*(1440) and 40-90 % for the A^*(1710)), 
it seems that they couple strongly to two meson-one baryon systems. There are many findings 
which support this idea, e.g., a strong aN coupling to the Roper resonance reported in [21} [22], 
an important contribution from the two meson cloud to the masses of the SU(3) antidecuplet 
members found in [23J, and a good reproduction of the data on the Svr distribution in the 
7r~p — > K^TjTt reaction by taking the ttttN decay channel of the A^*(1710) into account |24j . 
Hence, a study of the three-body structure of these resonances could shed more light on their 
properties. 

We study the TrvriV system by solving Faddeev equations in s-wave using the formalism 
developed in }25l [26l I27j . In these works, we studied systems made of two mesons and a 
baryon and those of three mesons by using unitary chiral dynamics to calculate the two-body 
amplitudes required for the Faddeev equations. In [25[ [26] we have investigated the ttKN 
system and coupled channels and found a strong coupling of the lowest lying 1/2"*" S and A 
resonances of the PDG [Ij to the three-body decay channels. In addition to this, unknown 
quantum numbers of some of the 5 = — 1 resonances were predicted, e.g., S(1560) has been 
listed with an unknown spin-parity in [1] and our work [25] generates it with J"^ = 1/2"*". 
Analogously, a three meson system formed by two pseudoscalar mesons, K and AT, and one 
vector meson, (/>, has been studied in [27], revealing the existence of a resonant state in this 
system at ~ 2150 MeV when the KK invariant mass is close to that of the /o(980), which 
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Figure 1: An example of a simplest possible interaction amongst the three particles, n^n^n. The 
labels ki (k'J on the particle lines denote the momenta corresponding to the initial (final) state. 
The meaning of the blob is shown in Fig. [H 



Figure 2: The blob in the Fig. \1\ which is a t-matrix. 



can be identified with the X(2175) 1 resonance discovered in the process e^e — > i;A/o(980) 
at BABAR [281 129] and in the J/* -^ ??(/)/o(980) reaction by the BES collaboration [30]. 

2 Formalism 

In this section we shall discuss different three body interaction diagrams, which contribute 
to the three-body scattering matrix, starting from the lowest order and going up to the 
higher order ones. We consider the -ktiN system with total charge zero and use vr'^vr'^n, 
vr+vr^n, vr^vr^n, 7r^7r~p and vr^vr'^p as coupled channels. Taking advantage of the formalism 
developed in [251 [26] , where 22 coupled channels were handled simultaneously, we have also 
used an extended base of states, up to 14, including vri^S, ttKA, TrrjN channels and we found 
the results remarkably similar to those obtained using the ttttN channels alone. Thus we 
present the formalism for the vrvrA^ states which is considerably simpler. We need the ttN 
and TTvr t-matrices as input, for which we solve coupled channel Bethe-Salpeter equations as 
explained in the following sub-sections. 

2.1 Lowest order diagrams 

Our aim is to calculate the three-body scattering matrix which includes all the possible 
"connected" diagrams between the three particles. The simplest possible diagram of this 
kind is the one which contains two i-matrices, for example the one shown in Fig. [T] for the 
ir^TT^n channel. Following [25|, this diagram can be expressed mathematically as ( reading 



Fig. [J from right to left as a convention 



t'g'h' = t' 
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where the superscript on t denotes the particle which is not interacting in the three-body 
system. Hence, t^ is the i-matrix for the interaction of particles 2 and 3 and t^ is that 
for particles 1 and 3. In order to calculate these vr^n — > 'K^n t-matrices the Bethe-Salpeter 
equation, 

t = v + vgt (2) 

is solved in the coupled channel approach with the potentials obtained from the Lagrangian 

[an [321 [Ml IS! 



Cmb = ^{Bij^'K^d^^ - d^m))B - B{^d^^ - d^^^))) 



(3) 



where / is the pion decay constant and the symbol ( ) denotes the trace in the flavor space 
of the SU{3) matrices $ and B 



$ 



B 



/ 1 0^ 1 

V2 Ve 


7V+ 

1 1 


K^ \ 


vr^ 


In 1 

V2 y/e' 




K- 
\ 


i?o 


"v/eV 


( >+ u 


V2 V6 


\ 

p 


s- 


n 


[I] 


E'^ 



(4) 



(5) 



Following [35j, vrA^, rjN, KA and KT, are taken as the coupled channels for the pion- 
nucleon system. For example, for total charge zero Eq. ([2|) is solved with the potential, 
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Figure 3: Different possible diagrams including two successive interactions. 



and the t^o„^^On element of the resulting matrix is used in Eq. ([T]) as t^ and t^. The two 
body propagator, g in Eq. ([2]), is divergent and is calculated using dimensional regularization 
by taking the substraction constants from [35], where the authors find the A^*(1535) as a 
dynamically generated resonance in the vriV system and its coupled channels. 

There are six possible three-body diagrams involving two t-matrices as shown in Fig. [3l 
To calculate all these diagrams, we require the tttt i-matrices also, which have been obtained 
by solving the Bethe-Salpeter equation (Eq. ([2])) with vrvr, vrr/ and KK as coupled channels 
[361 I37j . The potentials for these channels have been calculated using the chiral Lagrangian 
[m Ea [331 [36] 

Cmm = ^((5^$$ - ^d^^f + M$4), (6) 

ml 

M = \ ml 

2m|. - ml 

and niTT, m,^ are the pion and kaon masses, respectively. The two body propagator, g, in 
this case has also been calculated using the dimensional regularization [37]. A detailed study 
of these systems has been carried out in |361 [37] which revealed the dynamical generation of 
the a and /o resonances in the isospin zero sector and that of the ao in the isospin 1 sector 
of these mesons. 

All the diagrams in Fig. [3] can be expressed mathematically as fg^H^ with i 7^ j = 1, 2, 3. 
In the above discussion we have taken one channel, 7r'^7r''n, as an example but the calculations 
have been carried out by taking five coupled channels into account. Hence i*, g*-' and t^ are 



where 



matrices and each element of the g*-' matrix is given by 

S''ik^,k,)=(ll^]— -^, ^-^r-^ —l^i,l^j,= 1,2,3 (7) 



where D is the number of particles propagating between two t-matrices. Following the nor- 
malization of [38j, Nr = 1 for a meson and Nr = 2Mr for a baryon with Mj. being the mass 
of the baryon and ki (kj) is the momentum of the ith (jth) particle in the final (initial) state. 

2.2 Kinematics 

We now define the kinematics for the system. There are two variables in the calculation; the 
total energy of the three-body system, denoted as ^/s, and the invariant mass of the particles 
2 and 3, denoted as ^/s23■ The other invariant masses are obtained in terms of these variables 
as 

2 v"vv" ^ijy^'ZA I ""k ""jy 

.523 

with ?n,/j being the mass of the non-interacting particle and 



^ij = s + mi, (8j 



The definition in Eq. ([8]) implies an angular average between external momenta suited for 
the study of s-waves. 

From all this we can calculate the momenta, | fei | , | A;^ | of the particle 1 in the global center 
of mass system and that of the particles 2 and 3 in their rest frame (-R23)) which we denote 
as K {K') in the initial (final) state; 

|fci| = |Pi| = -^Ai/2(s,s23,m?) (10) 

2yJS 

\K\ = \K'\ = ^^X^I^'-- -2 ^2 



.. ,523,^2, m3J. 

^V'^23 

The calculation of the g^^ propagators for different diagrams requires the momenta of the 
particles in the global center of mass. For this, we boost the momentum in i223 to the global 



center of mass using the relations [39] 
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We define k\ to be along the z-axis and K to form a plane with ki, i.e., 








K 



\K\sin{6K) 




The f-matrices are calculated as a function of the invariant mass of the interacting parti- 
cles, e.g., t^ is calculated as a function of ^/s2z■ This is so because, in the chiral approach, the 



t-matrices can be split into an off-shell part, which behaves as {qj 



rrii 



with Qi being the 



four vector of the off-shell particle, and an on-shell part, where g^ is set to mf. In analogy to 
the findings in [25j, this off-shell dependence of the f-matrices in the three-body diagrams is 
found to cancel exactly with the three-body forces generated from the chiral Lagrangian for 
2 meson + baryon -^ 2 meson + baryon contact term in the SU(2) limit (see the Appendix). 
In [27j, where a study of the (j)KK system has been carried out and dynamical generation 
of the X(2175) resonance is found, an explicit calculation including the off-shell parts of the 
t-matrices has been done. The results obtained in this case are found to be qualitatively sim- 
ilar to those obtained by implementing the cancellation of the off-shell parts of the t-matrices 
with the chiral three body forces. The only difference in the results was a shift in the peak 
position in the squared amplitude by 40 MeV (about 2 % of the mass of the resonance). 
However it should be noted that the i-matrices and hence their off-shell parts, which have 
been derived from a chiral Lagrangian, are representation dependent. In the present work 
we rely upon the cancellation discussed above and retain the on-shell parts of the i-matrices, 
which depend on the invariant mass of the interacting pair in s-wave (as considered here). 

2.3 Higher order diagrams 

The calculation of the diagrams with more than two f-matrices involve a loop of three prop- 
agators and three two-body t-matrices. Such diagrams can be written as t^g'^H^ g^h\ where 
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Figure 4: A diagram involving three t-matrices. 



the propagators and the t-' -matrix depend on the loop variable even though the t-matrices in 
our model are calculated as a function of the invariant mass of the interacting pair. The cal- 
culation of the three-body scattering matrix would simplify if we could extract t^ g^H^ out of 
the loop integral. However, we would like to keep the loop dependence of this term. This can 
be done if the P g^^ calculated with off-shell variables is absorbed in the previous propagator 
as in [25j. Following the formalism developed in [25], we write the diagram shown in Fig. |4] 
as 

tlG^2l^2g2l^l ^ t\./J^)G''H\^J^)g^\k'„h)tHV^), (12) 

where 
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(14) 



Note that, while S23 is defined from the external variables for the diagram shown in Fig. 
m the argument S13 of the t^-matrix is a function of the loop variable and must be kept in 
the loop integral. We, thus, introduce t^ calculated as a function of 
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in F^^^ and hence in the loop integral G^"^^ . F^^^ also contains the inverse of t^ calculated 
as a function of S13 evaluated in terms of on-shell variables (Eq. dS])) and the g^^ propagator 
depending on off-shell variables along with the inverse of its on-shell version. 

In this way, [g"^"^ {k'2,ki)Y^ x [^^(a/sis)]""'^ in F^"^"^ (and hence in G^^^) when multiplied to 
t^{\/si^)g'^^{k'2, ki) in Eq. ()12p give an identity leaving t^g"^^ evaluated with the loop variable 
in the loop integral. Simplifying Eq. (|12p we have 
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Figure 5: A diagram with two concatenated loops. 



which is the right contribution of the diagram in Fig. HI The integrals of G*-' are regularized 
with a cut-off of 1 GeV in the modulus of the momentum, which if changed to 1.5 Gev 
introduces less than 1% of a change in G. One of the propagators in the G^^^ function is 
evaluated in the center of mass frame of two particles for convenience. A diagram with three 
t-matrices is, thus, written as t^G'^^ t^ g^ t instead of Vg^H^ g^ t . As has been discussed in 
[25], for diagrams with three t-matrices in general and in [26| for a specific diagram in detail, 
this scheme simplifies the calculations. The formalism is further developed by repeating the 
above procedure for higher order diagrams too, i.e., by replacing the g*-' propagator by the 
G^^ loop function everytime a new interaction is added. 

In case of the diagrams with more than three t-matrices, this scheme involves an approx- 
imation, since the introduction of a new interaction, to the diagrams of the kind shown in 
Fig. m replaces the external variables in the latter case by variables of a former loop. This 
procedure, which renders the integral Faddeev equations into a set of algebraic equations, is 
certainly very economical in terms of numerical solution and a justification for its use is given 
below. 

Let us discuss in detail a diagram with four t-matrices as shown in Fig. [5l as an example. 
This diagram is written explicitly as 
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which can be simphfied to 
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where s^^ is calculated analogously to Eq. (fT5]) . We compare our expression ([T8|) with the 
corresponding one written in terms of the g propagators of the concatenated two loops 
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(27r)3y (27r)3- 
The dependence of g^^ on the two loop variables has been thus factorized in Eq. (jlSp as 

g^Hquq2) = ^1(51)^2(92), (20) 

where 



and 



=^2(92) =5 (9*2 
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(22) 



2Ei{q2) E?Xq2) Vsi3 - ^1(92) - -^3(9*2) 

This factorization, which simplifies the calculations to a great extent, leads to very similar 
results to those obtained with the concatenated loop function as can be seen in Fig. [6l where 
we show the mod-square of the Eq. (llSp and Eq. ()19p as a function of \/i, in the energy 
region of our interest. The agreement of the results depicted in Fig. [6] shows that Eq. (jlSp 
is a good approximation of the Eq. (I19p . Hence, this scheme is used to write the rest of the 
higher order diagrams which contribute to the three-body amplitude. 

If we sum Eqs. (P), (fT2]) . (fT7|) and all the other possible diagrams with the last two 
t-matrices as t^ and t^, we get the series 



t'g'H' + t'G''H'g'H' + eG''H'G''H'g'H' + ■ ■ ■ + t'G''H'g'H'' + t'G''-H'G'^H^g^H' + • • • , 

(23) 
which we define as Tj^. Similarly, we consider all other possible diagrams obtained by per- 
mutating different interactions between the three hadrons and get the following equations 
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Figure 6: The comparison of modulus square of Eq. U^) and Eq. [W\) shown as dashed and solid 
lines, respectively, in units o/10~^^ MeV~^ . 



upon suniniing all the diagrams with the same last two t-matrices 
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These are six coupled equations which are summed to get 

Tn = tA^ + rA=^ + rf + rf + rf + rf . 

The Tjj can be related to the Faddeev partitions T* of the the Faddeev equations 
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3 Results and discussions 

The Tr in Eq. (j25p has been calculated in s-wave for the coupled channels TT^TT^n, n^-K^n, 
7r~7r~^n, n^n'p, tt'tt^p as a function of ^/s and y^sii- All the angle dependent expressions 
have thus been projected into s-wave. The Tfj-matrix (Eq. (|25p ) is then projected on the 
isospin base defined in terms of the total isospin of the three body system, /, and the total 
isospin of two pions, 1-,^, defining the states as |I, /^tt)- These states are obtained assuming 
the phase convention for | vr"*") as — | 1, 1). We write the state | tt^ vr*^ n), for example, as 

I 7r° 7r%) = I 1, 0)(g) I 1, 0)0 I 1/2, -1/2) 

l\U^ = 2,1^^ = 0) - y^ I /.. = 0,1^, = 0) U 

I 1/2,-1/2) 

^ I / = 5/2,/^, = 2) + -^\I = 3/2,1,^ = 2) - y 1 I / = 1/2,1,^ = 0) 



To simplify the notation, we omit the label / and /tttt and write 

I ^Ovr^n) = y^ I 5/2,2) +-^ I 3/2,2) - ^^ | 1/2,0). (28) 



Similarly, 

k+vr-n) = -yi|5/2,2)-yi|3/2,2)-yi|3/2,l)-yi|l/2,l)-yi|l/2,0) 

Ivr-vr+n) = -^1 | 5/2, 2) - ^^ | 3/2, 2) + ^^ | 3/2, 1) + ^^ | 1/2, 1) - ^i | 1/2,0) 

|vr-7r%) = yi|5/2,2)-yA|3/2,2)-yi|3/2,l)+yi|l/2,l) 

|vr°vr-p) = yi|5/2,2)-y^|3/2,2) + yi|3/2,l)-y^|l/2,l). (29) 



From Eqs. (|29p . one can obtain 



5/2,2) = W-(\/2 I 7r°7r°n)+ I 7r°7r"p)+ I vr^vr^p) - W- I TT+vr^n) - W- I TT^vr+n) j 

/I / 3 3 

— 2 I 7r°7r°n) - ^ I 7r°7r"p) - ^ I tt'tt^p)- \ Tr+vr^n)- | tt" 7r+ n) 
15 \ \/2 v2 

1/2,0) = -W-M 7r°7r°n)+ I 7r+7r"n)+ I TT^vr+n) 
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Figure 7: The squared amplitude for the irnN system in isospin 1/2 configuration as a function of 
i/s and y/s^. 



One could equivalently define the states in terms of the total isospin and the isospin 
of a pion-nucleon sub-system {Ittn) by repeating the former procedure or using the Racah 
coefficients for the transformation of the | I, Inn) states to | I,InN) states. 
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In Fig. [7] we show the squared amplitude | T^ 



R 



Tr 



Y, t'g'Jp 1 2 for the vrvriV 

system, calculated in s-wave and projected on the isospin base | I, Inn) =| 1/2,0). The Y1^9^ 
has been subtracted out of the Tr (following [2^) since it does not give rise to any three-body 
structure and only provides a background to the amplitude. The squared amplitude shown 
in Fig. [7] has a peak at ^/s = 1704 MeV, with a full width at half maximum of 375 MeV (see 
also Fig. [8l). These results are in good agreement with the characteristics of the A/^*(1710) 
[1] and, hence, we relate the resonance shown in Fig. [7] with the A^*(1710). To get further 
physical meaning of this peak, we show the same amplitude depicted in Fig. \7\ but as a 
function of ■^/s23 and y^sU in Fig. [9l The peak in ^/s^ is very wide (width ~ 270 MeV ) and 
is in the energy region of the a resonance (see also Fig. [TUj) . This means that the A^*(1710) 
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Figure 8: The projection of the amplitude shown in Fig. [^ 
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Figure 9: The same as shown in Fig. \7\ but as a function of the tttt invariant mass and that of the 
TiN system. 
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Figure 10: The projection of the amplitude shown in Fig. 



has a large vrvrA'^ component where the tttt sub-system rearranges itself as the a resonance. 

Although we find evidence for the A^*(1710), this work fails to find any clear trace of the 
Roper resonance, which means that considering the vrvrA^ system in s-wave interaction does 
not suffice to generate the Roper resonance, which is not surprising. Other works such as 
the Juelich model [21], which successfully describes the dynamical generation of the Roper 
resonance, contains additional information on the irN, TrA,pN coupled channels and aN 
forces beyond the three body contact term of the chiral Lagrangians which we include here 
and which cancels the off-shell dependence of the amplitudes. An important contribution of 
the vrA channel and vrvr final state interaction (with one of the pions coming from the decay 
of the A resonance) to the Roper resonance has also been claimed in [40j . Such information 
is not present in our formalism. Things are different in the case of the A^*(1710) with its 
large empirical coupling to vrvrA^ and weaker to vrA^ and other coupled channels. 

Another important result of this work is that we do not find any resonant structure in the 
total isospin / = 3/2 and / = 5/2 configuration. Should we have found the latter, it would 
be exotic in the sense that it would not be possible to construct it with just three quarks. 
But no structure is found in this isospin state. 

To summarize, we have studied the vrvrA^ system in s-wave, thus in J^ = 1/2+ configu- 
ration. We find a resonance, in three-dimensional plots of the squared amplitude versus the 
total energy and the invariant mass of a sub-system, at 1704 MeV, which can be associated 
with the iV*(1710) [ij. Our peak has a full width T = 375 MeV to be compared with that of 
the A^*(1710) which ranges from 90-500 MeV [1]. We find that the invariant mass of the vrvr 



15 



sub-system falls in the region of the mass of the a (500 -i 200 MeV) when the vrvrA^ amplitude 
peaks at ^/s = 1704 MeV, which means that the large width of the A^*(1710) could be related 
to that of the a resonance formed in the vrvr sub-system. No evidence for states with / = 3/2 
and / = 5/2 is found in this work. We also do not find the Roper resonance in our approach. 
This should not be seen as a negative result, but as an evidence that the structure of the 
Roper is far more complex than that envisaged by the vrvrA^ interaction in s-wave, which is 
what we have investigated in the present work. 
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Appendix 

In this appendix we discuss the three-body interactions including the off-shell parts of the 
t-matrices that give rise to a kind of three body force which, as we show below, gets cancelled 
with the three-body force arising from the chiral Lagrangian. 

Let us consider the lowest order diagrams, which correspond to the first terms of the Tr 
equations, i.e., t^g^H^ . There are six terms of this kind shown in the Fig. [3l which can be 
expanded in terms of the potentials as 



fg'-H^ = [v^ + v'-g'-v'- + v'-g^v^g'-v^ -\ ] g"-^ [v^ + v^g^v^ + v^gH^g^v^ -\ ] 

= v'^-g'-^v^ + v'-g'v'g''^v^ + v"- g'-' v^ g^ v^ H (30) 

For example, a term of t^g^'^t^ expanded as in Eq. (j30p is shown in Fig. [TTJ 

The potentials in chiral dynamics can be split into an on-shell part which depends on 
the center of mass energy of the interacting particles and an off-shell part proportional to 
p'^ — m^ for each of the meson legs, in case of meson-meson interaction (where p is the four 
vector of the off-shell particle and m is its mass). In case of the meson-baryon interaction, 
the off-shell part of the potential behaves as p^ — k^ , where p^{k^) is the energy corresponding 
to the off-shell (on-shell) momentum. Due to this behavior, the off-shell part of the potential 
cancels a propagator in the loops, giving rise effectively to a three body force, for example, 
the one shown in Fig. [T2] corresponding to the t^g^^t^ term shown in Fig. [TTl 

Similar effective three-body forces arise from other terms too. We shall now write the 
contributions for the first terms of all six fg^H^ terms (Eq. ()30p ) including the off-shell parts 
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(a) 



(b) 



Figure 11: A diagrammatic representation of the t^g^^t^ term. The blob in (a) represents a t-matrix 
which can be expressed mathematically as v + vgv + vgvgv + vgvgvgv + .... And (b) shows the term 
{v^g^v^g^v^)g^^{v'^g^v^g^v'^) of Eq. (E^)- 



Figure 12: An induced effective three-body force generated by the cancellation of the off-shell part 
of the potential and a propagator as explained in the text. 

of the t-matrices, taking the vr^vr^n channel as an example and evaluate the total effect of 
these three-body forces. 

We label the initial (final) four-momentum of the vr"*" as p {p'), that of the tt~ as k {k') 
and that of the neutron as q {q') as shown in Fig. [13l We assign a four vector k" to the 
intermediate states, see Fig. [TT] 

7r+ {p) TX+ {p') 

TT" {k) \ \ -K' [k') 

\ ] 

n(g) Vv n{(^) 

Figure 13: Assigning four momenta to the tt^tt^u system. 

The potentials calculated from the chiral Lagrangians Eqs. ([31 [6]) for the three possible 
two-body interactions are 



K+ 



1 



3s^^ - ^{Pi - mf) , (31) 



1 '-.0 , J./0^ 



Kr-n->7r-n — ~T72 ^ t + ^tt ); (32) 

17 



K 



7rf n^TrT"?! 



-K 



TV n — ^■Tr Wi 



(33) 



where Stttt is the invariant mass of the vr — tt sub-system, / is the pion decay constant and 
k^ {k'^ ) is the energy of the pion before (after) the vrA^ interaction. 

In this way, the contribution of the first term of Eq. (j30p for i = 1 and j = 2, which 
corresponds to the first diagram in Fig. [3l is given by 



' {k' + k"^"" 



1 



16/ 



Enk'^ + q'^ - fcO - Enip' + k) + 



-(pO+p'°) = T- 



(34) 



le 



rUn in Eq. (j34p is the neutron mass and the superscript "on" on Tq denotes that there is no 
off-sheh dependence in the above equation. 
For the diagram (b) of Fig. [3] 



n 



1 



,//0 



u/On 



24/4 L 



2k"' + {k"" - k'") 



1 



k" — mt 



2,{p + kf 



t -ml 



rpon ^ j^off 



with 



rpon 

T°ff 



':k" 



4/4 [p + k-p 



/^2 



mt, 



-{p + kf 



24/4 



k 



/O 



-p^ -k^ +p'^ + ^{p + k) 



,k"^-k'' 
k" — m'i 



(35) 
(36) 

(37) 
(38) 



representing the on-sheh and off-shell contributions to Th. In Eq. (j37p and in the first term 
of Eq. (|38p . k" has been replaced by p + k — p' using the energy- momentum conservation 
law from the initial state. For the second term of Eq. (j38p we apply energy-momentum 
conservation from the final state 



. //2 



k"^ = (k' + q'-qf. 
Defining Aq = q' — q, Eq. ([39]) becomes 



(39) 



and, hence 



Therefore, 



k 



k"=ml + {Aqf + 2k'-Aq 
"0 f^/0 f^^_-^o 



(Aq)' 



k"^ 



mt, 



b 24/4 



A:'°-pO 



(A^)2 + 2A;'- A^' 



,/0 



k" + p"" + 3{p + ky 



(Ag)0 



(Ag)2 + 2k' ■ Aq 



(40) 



(41) 



The contribution of the diagram (c) of Fig. [3] is 



16/ 



1 / , /0\1^n 
— (P +P ) 



uO , , /O 



Enk^ + qO-k'^-Enip + k') 



(A;V/fc')=Tr. (42) 



Similarly, 



Trf 



1 



/O , //O /Ox 



1 



(2p +p -p' J 

p" — m^ 

on I rroff 



2 ^ //2 _2n 



24/ 



[3(p + A;)^-(/^-m^)] 



with 



rpon 



rpOff 
-'d 



1 /O 



4/4^ (p + A;-A:')2-m: 



(p + A:)^ 



24/4 



p' — m^ 



Analogous to Eq. ()40p . we write 



P"'-P" 

V'2 _ TTi2 



(Ag)C 



y/2 _ ^2 (Ag)2 + 2p' • Aq 



which gives 



t: 



'Off 



1 



'0 „0 7.0 , , /O 



24/4 
For the next diagram, we have 



p'-^p^-k^ + k'^' + Sip + kf 



(Ag)0 



(Ag)2 + 2p' • Aq 



^ [3{p' + k'f - ik'" - ml)]-^^{2k' + k"' - k') 



24/ 

rpOn _|_ rpOff 



' k" — m^ 



where 



rpon 
J'off 



' (p' + '^rj-T-T, 



4/4 
1 



-A;° 



(p' + /c' — p)2 — ?n,2 
A;0_p'0_A;'°+p° + 3(p' + A: 



,0 „/0 7/0 , n , o/„/ I 7 /\2 ^ ^ 



fc"^ — m^ 



In this case k" = k — Aq and therefore 



(Ag)C 



fc'/2 _ ^2 (Ag)2 - 2A: • Ag 



(43) 
(44) 

(45) 
(46) 



(47) 
(48) 

(49) 
(50) 

(51) 



(52) 
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Figure 14: Source of three-body force from the chiral Lagrangians. 



leading to 



-o// 



24/^ 



uO „/0 , /O , n 



k^_p'^ _k"'+p^-?,{p' + k') 



i\2 



(Ag)0 



(Ag)2 -2k- Aq 



For the last diagram of Fig. [3] we have 

1 



f/ 



24/^ p" — m^ 



(2p0+/°-p0) 



rpon _, rpOff 



where 



rpon 



rpOff 



4/4^^ (p' + /fc'-/fc)2-m2 



P 



1 



24/4 
Following the same method 



.0 ^/O , /O I ,0 



pO-p'"-A:'VA:0 + 3(p' + A;') 



p" — m| 



(A^)C 



then 






24/4 



(Ag)2-2p- Ag' 
pO - p'O - k'"" + k"^ - 3{p' + k'f 



{Aqf 



(Ag)2 -2p-Aq 



(53) 

(54) 
(55) 





(56) 




(57) 



(58) 



(59) 



On the other hand, genuine three-body forces also originate directly from the chiral La- 
grangian , where we can find a contact term as the one shown in Fig. ()14p |41j. 

At lowest order in momentum, which we consider in our study, the interaction Lagrangian 
between mesons and baryon is given by 



C = i{B^>^[T^,B]) 



(60) 
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where 



T^ = \{u^d^u + ud^u^), u^ = e'^^^'i 



(61) 



and 0, B are same as those in Eqs. ([H [5]). If we expand T^ up to the terms which contain 
four meson fields, we get 



r„ 



1 






For the case under consideration, i.e., 7r"'"7r~n, the Eq. (j60p becomes 



(62) 



£ 



32/' 



-n 



—^ir TT^vr vr^ — vr gf-K^ir vr^ + vr vr^^vr vr^ vr vr^vr ^vr^ 



n (63) 



(64) 



In this way, the contribution of the diagram in Fig. [T3] is 

Tsb = ^Mq'X'^^- 2^ - 2^+ 2/)n,(gl 

We are interested in the low energy region, thus, only the 7'^ component of Eq. (j64p is 
relevant, then 



1 



f36 



24/' 



-(2pO-2A;'°-2A;0 + 2p'°) 



(65) 



Adding this to the off-shell contributions from the Faddeev equations at second order in 
t-matrices, we get 



E^°^^+^3. 



1 



i=l 



24/4^ 

+3{p + k)\Aqf 



4fc° + 4p'°-4A:'° + 4/ 
1 



1 



+3{p' + k'y{Aq) 



(Ag)2 + 2/c' • Ag {Aq)^ + 2p'-Aq 
f 1 1 



(Ag)2 -2p-Aq (Ag)2 -2k- Aq 



(66) 



If we consider small momentum transfer for the baryon, i.e., A^ << 1, Eq. (j66p can be 
expressed as 



Et-" + t,, 



1 



4 = 1 



24/4^ 
+3ip' + k'f 



4A;° + 4p'° - 4A:'° + 4p° + 3{p + kf 
' 1 1 ' 



1 



1 



(Ag)0 - 2p0 (Ag)0 - 2A;0 



(Ag)0 + 2A:'0 {Aq)0 + 2p'^j 

(67) 
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And there is a cancellation of the terms in the SU{2) limit, assuming equal average energies 
for the pion. Furthermore, if the propagators in the Eq. (|66p are projected over s-wave, as 
we do in our study, the curly brackets become 

f 1 , f (Aqf + 2k'\Aq)^ + 2 \ k' \\ Aq 

-^ — ln\ 



2\k' \\Aq\ \{AqY + 2k'^{Aqf -2 \ k' \\ Aq 

1 , ({Aqf + 2p'\Aqf + 2\p'\\Aq 

-^ — ln\ 



and 



2\p'\\Aq\ \ (A^)2 + 2p'^{Aqf -2\p'\\Aq 



1 , / (Aqf - 2p^(Aqf + 2\p\\Aq 

-Inl 



2\p\\ Aq\ \ {Aqy - 2p^{Aqf -2\p\\Aq 

1 , ( (Aqf -2k^iAqf + 2\k\\Aq 

-ln\ 



2 I A: I I Ag I \ {AqY - 2k^{Aqf - 2 | A; | | Ag 

respectively, and the cancellation is exact. 

With the cancellation of these basic diagrams proved, the addition of an extra interaction, 
t;*, to the set of these cancelling terms will still give a vanishing contribution. For the on-shell 
part of the basic diagrams that we have studied, further iteration of the potential between 
two particles, leading to the two-body t-matrix, are done in such a way that the off-shell part 
of the new v"^ interactions is reabsorbed in constants of the on-shell potential, as is done in 
the construction of the two body i-matrices [Ml ES] • The conclusion is that at the end we 
should only use the on-shell t-matrices, ignoring the off-shell effects and genuine three body 
forces simultaneously. 
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